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Gerard  E.  Welch*  and  Louis  M.  Larosiliere* 

Army  Research  Laboratory 
NASA  Lewis  Research  Center 
21000  Brookpark  Road,  M/S  77-7 
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The  unsteady  flow  within  wave  rotor  passages  is  influenced  by  the  rotor  blade,  hub,  and  tip-shroud 
surface  profiles.  By  averaging  from  hub  to  shroud  and  from  blade  to  blade,  a  reduced  set  of  the 
governing  equations  is  obtained  that  is  appropriate  for  design  studies  and  parametric  analyses.  The 
application  of  these  equations  requires  closure  models  for  force  integrals  and  for  correlation  terms  that 
arise  when  the  density-averages  of  products  of  the  flow  field  variables  are  expanded  in  terms  of 
products  of  the  density-averaged  variables.  The  force  integrals  and  the  correlation  terms  depend  on 
the  instantaneous  pitchwise  and  spanwise  flow  field  distributions  established  by  unsteadiness  relative 
to  the  rotor,  flow  turning  induced  by  blade,  hub,  and  tip-shroud  profiling,  and  rotation.  Two 
approaches  to  model  the  force  integrals  are  described.  The  influence  of  relative  unsteadiness  and  flow 
turning  on  the  correlation  terms  is  discussed  by  considering  the  propagation  of  gas  dynamic  waves  in 
rotor  passages  defined  by  uncambered,  staggered  blades  and  by  unstaggered,  cambered  blades. 


Nomenclature 

1/2 

=  iyplp)  5  speed  of  sound 
e  -  unit  vector 

=  -  (p/  p ) ,  specific  relative  energy 

=  inviscid  flux  vector 
=  flux  vector  of  correlation  terms 
=  source  vector 

=  correlation  terms  of  radial  source  terms 
=  A  ~  r  Q  M  ,  specific  rothalpy 
=  — -pip  +  -(w-w),  specific  total  enthalpy 

=  vector  of  force  integrals 
=  axial  length  of  rotor 

=  rotor  chord  length 

=  rotor  pitch 

=  static  pressure 

=  vector  of  conserved  variables 
=  (Ry[i),  polar  coordinate  system  of  Fig.  3 
=  (r,0,z),  position  vector  relative  to  spinning 
rotor 

S  =  passage  cross-sectional  area  (cf.  Eq.  5). 
t  =  time 

u  =  absolute  velocity  (e.g.,  )) 

w  =  M  -  Qre^,  relative  velocity  (e.g., 

w  =  relative  speed  of  shock  wave  (see  fig  1). 


z*  -  normalized  axial  distance  (see  Figs.  1-3) 

P  =  tan  (w^/w  ),  relative  swirl  angle 

Y  =  ratio  of  specific  heats 

A  =  A  ((p)  =  <p  -  <p  ,  pitchwise  difference 

A  =  A((p)  =  (p-(p,  spanwise  difference 

e  ^  =  rafio  of  static  pressures  across  expansion  fan 
o  =  solidity 

5  S  j  “H )  9  intrinsic  coordinates 

p  =  mass  density 

4)  =  tan"  (w  /w  ),  meridional  flow  angle 

Q  =  shaft  angular  speed 

=  denotes  an  area-average  (cf  Eq.  5) 

<>  =  denotes  a  density-weighted  average  (cf.  Eq.  6) 

Subscripts 
H  =  rotor  hub 

P  =  rotor  passage  leading  blade-surface 

S  =  rotor  passage  trailing  blade-surface 

T  =  rotor  tip-shroud 

0  =  rotor  passage  initial  relative  conditions 

1  =  rotor  passage  inlet  (or  exit)  relative  conditions 

Introduction 

The  wave  rotor  is  a  member  of  the  family  of 
dynamic  pressure  exchange  devices.^  It  consists  of  a 
shrouded  rotor  surrounded  by  a  stationary  casing.  The 
casing  endwalls  are  penetrated  in  sectors  by  inlet  and 
outlet  ducts  that  port  gas  to  and  from  the  rotor  flow 
annuli.  Gas  dynamic  waves  are  initiated  and  terminated 
as  the  rotor  passages  open  and  close  to  the  ported  flows. 
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Wave  rotors  are  self-cooling  and  can  therefore  be  used 
to  top  turbine-inlet-temperature-limited  gas  turbine 
engines.^  Recent  research  has  emphasized  axially 
bladed  wave  rotors  (Le.,  pressure-exchangers);  however, 
engine  performance  levels  might  be  enhanced  more 
significantly  by  topping  with  wave  rotors  that  produce 
net  shaft  power  (Le.,  wave  engines  or  wave  turbines). 

Motivated  by  this  potential  advantage,  a 
design/analysis  tool  is  currently  needed  which  is  able  to 
predict  accurately  the  performance  of  wave  rotors  with 
arbitrary  passage  geometries.  One-dimensional  wave 
rotor  flow  solvers  suitable  for  parametric  analyses  can 
be  used  to  predict  the  performance  of  pressure- 
exchangers  accurately;^  however,  there  are  no 
comparable  computational  tools  available  for  the 
generalized  rotor  passage  geometries  of  wave  engines. 
Although  multi-dimensional  flow  solvers^*'^  can  be  used 
to  compute  the  wave  rotor  flow  field  in  complicated 
geometries,  it  is  currently  impractical  to  perform  the 
parametric  analyses  of  design  studies  with  these  tools. 

The  purpose  of  the  present  work  is  to  formulate  a 
passage-averaged  description  of  wave  rotor  flow  which 
accounts  for  blade  stagger  and  camber,  and  hub  and  tip- 
shroud  profiling  in  a  reduced-dimension  model  suitable 
for  parametric  analysis  tools.  The  paper  proceeds  as 
follows:  1.)  the  passage-averaged  formulation  is  obtained 
by  pitchwise-  and  spanwise-averaging  of  the  Euler 
equations;  2.)  closure  models  for  the  blade  force 
integrals  resulting  from  the  averaging  procedure  are 
provided;  and  3.)  the  influences  of  relative  unsteadiness 
and  relative  flow  turning  on  correlation  terms  introduced 
in  application  of  the  passage-averaged  equations  are 
discussed  by  using  analysis  of  the  flow  fields  established 
by  the  propagation  of  gas  dynamic  waves  in  idealized 
wave  rotor  passages. 

Passage-Averaged  Model 

Wave  rotors  are  currently  envisaged  to  be  axial  or 
mixed-flow  machines.  Influenced  by  the  rotor  geometry, 
the  Euler  equations  are  averaged  fi*om  hub  to  tip  and 
from  blade  to  blade  in  this  work.  This  averaging 
procedure  yields  a  set  of  equations  with  one  (axial) 
space  derivative  suitable  for  a  CFD-based  wave  rotor 
analysis  tool. 

Governing  Equations 

The  mass,  momentum  (component  in  an  arbitrary 
direction  e.) ,  and  energy  equations  for  an  in  viscid  flow 
written  in  a  frame  of  reference  that  is  rotating  with  the 
rotor  are 

—  +  V-(p\v)  =  0  (1) 


3p(w*e.) 

- +V-(pw(w*ep) 

Di, 

=  p[(Q^r  -  2Q ^w)-e.  +  w*-^] 


(2) 


and 

dpi 

- -  +V-{pwfi  )  =  0  ,  (3) 

respectively.  In  cylindrical  coordinates,  for  example, 
De  IDt  =  w  e  Ir  ,  DeJDt=  -  wi  /r,and  Di  IDt  =  0. 

T  0  0  8  0  r  2 


Averaged  Equations 

By  averaging  Eq.  1  from  hub  to  tip  and  from  blade 
to  blade  in  cylindrical  coordinates  it  is  found  that 


dps  5pS<w^> 
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dt 


dz 
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(4) 


where  the  area-  and  density-averages  are  defined  by 
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and  where  flow  tangency  has  been  imposed  at  the  hub, 
tip-shroud,  and  blade  surfaces.  S  here  denotes  the 
passage  cross-sectional  area. 

By  applying  the  averaging  procedure  to  the 
momentum  equations  (Eq.  2)  in  the  relative  cylindrical 
coordinates,  the  following  equations  are  obtained: 


(Axial-momentum) 


a 

di 


(p5<w>)  +  — (p5<ww>+p5)  = 

z  dz  2  ^ 

0  ^ 

P  T 

j A  (prtan(j))d0  +  J A^(pim^)dr 


(7) 
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(Tangential-momentum) 

a  _  3  _ 

—  (pS<ur>)+  —  {pS<wuj>)  = 
dt  9  dz 


Z  0 


(8) 


By  using  these  relationships  to  expand  the  density- 
averaged  products  in  the  conservation  equations  above, 
the  passage-averaged  equations  are  obtained: 

dQ  dF  dF'  , 

^  +  _z:  =  _=  +  G  +  G  *L 
dt  dz  dz  ~  ~  ~ 

where 


(15) 


(Radial-momentum) 

3  _ 


-  -  3  _ 

—  (pS<w>)  +  —  (p5<ww>)  = 
dt  ’■  dz  ^ 
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-  JA  (j>r)dB  +  p5  < 


pip 


(9) 


where  tarKj)  =  dr/dz  on  the  hub  and  shroud  surface  and 
tanp  =rdQldz  on  the  blade  surfaces.  It  was  deemed 
appropriate — compare  the  time  derivative  and  the 
convective  flux  terms — to  express  the  flow  field 
variables  in  terms  of  density-weighted  averages  and 
area-averages.  Note  that  the  averaged  equations  have 
been  expressed  in  terms  of  the  absolute  velocity 
M  =  w  +  r  Q  e  and  the  momentum  equations  are 
influenced  weakly  by  the  rotor  speed  through  the 
equation  of  state 

<:^>  =  (Y-l)(<c^>-i[<w-w>-<(rQ)^>])^10) 
where  <pl  p>  =pl  P  ^ 

Application  of  the  passage-averaging  procedure  to 
the  relative  energy  equation  (Eq.  3)  provides 


dQS<e >  dpS<w  h  > 

- fj_  =  0 


(11) 


dt 


dz 


Relating  Averaged-Products  to  Products  of  Averages 
In  application  it  is  necessary  to  expand  the  averages 
of  variable-products  (e.g.,  <w  w  >)  into  products  of 
averaged- variables  (c.g.,  ><w  >).  A  deviation  from 

the  mean  of  an  arbitrary  vanable  /  is  defined  as 


P/ 


/./-</>  =/-^ 


(12) 


By  applying  the  passage  averaging  procedure,  the 
following  product  relationships  are  obteiined: 

(13) 
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and 
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The  blade  and  hub  and  tip-shroud  force  integrals 
must  be  modeled  in  order  to  apply  the  passage-averaged 
equations.  The  correlation  terms  too  should  be  modeled 
or  else  shown  to  be  negligible. 


Force  Integrals 

The  momentum  equations  are  weakly  coupled 
through  the  blade  and  hub  and  tip-shroud  force  integrals 
of  Eq.  21.  The  force  integrals  are  related  to  the  pressure 
gradients  sustained  by  acceleration  induced  by 
unsteadiness,  flow  turning,  and  rotational  forces.  Two 
approaches  to  close  the  system  of  equations  are 
suggested  here:  a.)  a  local  flow  field  can  be  postulated 
for  the  sake  of  deriving  expressions  for  the  force 
integrals  of  the  three  momentum  equations;  and  b.)  the 
three  momentum  equations  can  be  collapsed  through  the 
force  integrals  into  one  effectively  axial  momentum 
equation. 

Local  Flow  Sub-Model  Approach 

Consider  the  case  of  a  gas  dynamic  wave 
propagating  normal  to  the  blades  of  a  staggered, 
uncambered  passage  (see  Figs.  1  and  2).  For  the  sake 
of  brevity  in  this  discussion,  radial  variations  are 
neglected  and  only  the  blade-to-blade  (tangential)  sub¬ 
model  is  considered  here.  Because  r30/dz  =  O, 
dpldr\  =0.'  The  local  loading  is  therefore  related  to  the 
axial  gradient  of  the  area-averaged  pressure  without 
approximation  through 

dpS 

Ag(/»)  =  Pj, -P^  =  sinp  cosP  —  (22) 

where  P  is  used  here  to  indicate  the  constant  stagger 
angle.  The  axial  and  tangential  momentum  equations 
can  then  (in  this  very  restrictive  case)  be  rewritten  as 

3  _  3  — 

—  (pS<>v  >)  +  — (piS<w  w  >)  + 
dt  z  dz  ^  ^  (23) 

—  (p5cos*P)  =  0 
dz 


a  _  a  _ 

—  (pS<M  r>)  +  —(pS<w  u  r>)  + 

a,  >  az  .  •  p4) 

—  (p5rcosPsinP)  =  0  . 
dz 

Although  the  sub-model  approach  is  intended  to  retain 
the  various  momentum  equations,  note  that  if  it  is 
further  assumed  that  locally  =  w  tanP ,  then  Eq.  24 
is  redundant  and  simply  re-expresses  continuity  and 
axial  momentum  conservation. 

It  is  beyond  the  scope  of  the  present  paper  to 
develop  a  detailed  sub-model  for  the  force  integrals; 
however,  two  ideas  are  formed  from  the  above  results. 
First,  it  is  possible  (in  at  least  one  case)  to  develop  a 
sub-model  for  the  local  pressure  distribution — in  this 
example,  dp /dr]  =0 — which  can  be  used  to  relate  the 
local  loading  to  the  axial  gradient  of  area-averaged 
pressure.  (Future  work  would  be  required  to  develop  a 
generalized  local  flow  sub-model,  although  it  is  unclear 
whether  or  not  this  a  practical  goal.)  Second,  the 
system  of  equations  is  reduced  by  relating  the  tangential 
and  radial  velocity  components  to  the  axial  velocity 
component.  This  approach  is  used  below. 


*'One-Dimensionar*  Model  Approach 

By  relating  the  density-averaged  tangential  and 
radial  velocity  components  to  the  density-averaged  axial 
velocity  component  via  specified  area-averaged  flow- 
tangent  axial  distributions,  the  axial  loading  term  can 
be  related  to  the  tangential  and  radial  loading  terms 
and  L  ,  and  the  three  momentum  equations  can  be 
folded  into  a  single  axial  momentum  equation.  If  the 
correlation  terms  are  neglected  (for  the  moment),  then 
the  blade  force  integral  of  the  axial  momentum 
equations  can  be  expanded  into 


y* A^(/7tanp)dr  «  tan|J  jA^(p)dr  + 


(25) 


pj  Ag(tanP)dr 


and  the  hub  and  tip-shroud  force  integral  can  be 
expanded  into 


J  A  «  tan(^  J  A^(j?r)dd  + 

6 

p 

pr  J  A  (tanp)rf0 


(26) 


where  it  is  also  assumed  that  pr  -  p  r.  The  momentum 
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equations  are  related  through  the  force  integrals  and  can 
be  folded  together  by  assuming  that 
<w  > 

— —  =  tanp  “  -(tanPI  +tanP|  )  (27) 

<w  >  2  P  S 

z 

and 

<w  > 

— 1-  =  tan$  “  -(tancj)!  +tan4)|  )  .  (28) 

<>V  >  2  T  H 

z 


It  is  noted  that  the  error  introduced  in  these  assumptions 
is  likely  on  the  order  of  the  error  incurred  by  neglecting 
the  correlation  terms.  Under  these  assumptions,  the 
passage-averaged  equations  can  be  written  as 


pS 

pS<w> 

0  ‘ 

d 

pS<w  >6^ 

d 

pS<w  +  pS 

G 

(29) 

Jt 

z 

pS<e> 

dz 

Z 

p5<>v  ><h> 

z  i  \ 

0 

where 


6^  =  1  +  tanp^^  +  tan(|)^  . 


and  where 


“  p  5  [  tantj)  Q  <r>  2  dz 


tanS 

pS — I  +p5[— - 

<r>  r^-r 


A  (tan(j))  A  rtanp) 

r  0 


r(e  -0  ) 

T  H  ^  P  s' 


(30) 


]  + 


(31) 


Equations  29-31  represent  a  closed  set  of  effectively 
one-dimensional  equations  that  account  for  the  influence 
of  blade,  hub,  and  tip-shroud  profiling.  The  source 
term,  Eq.  31,  evidently  contains  the  influence  of  rotation 
and  blade,  hub,  and  tip-shroud  surface  profiling  on  the 
rotor  flow.  Note  that  by  applying  Eq.  29  in  the 
restricted  c^e  of  a  passage  with  constant-radius  hub  and 
tip-shroud  surfaces  and  uncambered  blades  at  stagger 
(i.e.,  6"  =  l/cos"p  and  G  =0),  Eq.  23  is  recovered. 
Further,  the  classical  quasi-one-dimensional  equations 
are  recovered  in  the  limiting  case  of  8  =  1  and 
tan(^  =tanp  =0. 

Correlation  Terms 

The  correlation  terms  (e.g.,  <w  w  >)  were 
introduced  by  the  expansion  of  the  density-average  of 
variable-products  (<w  w  >)  into  products  of  the 

density-averaged  variables  {<w  ><w  >).  These  passage- 
averaged  "stress"  terms  are  non-zero  if  there  are 
instantaneous  pitchwise  or  spanwise  variations  in  the 


product  constituents  (w  and  w).  The  variations  are 

induced  in  general  by  a.)  acceleration  caused  by  relative 
unsteadiness,  relative  flow  turning,  rotational  forces,  and 
shear  forces  and  b.)  entropy  nonuniformities  due  to  the 
presence  of  shock  waves  or  contact  discontinuities 
oriented  other  than  normal  to  the  z-axis,  heat  transfer, 
and  wakes.  Although  the  local  values  of  the  correlation 
terms  may  be  large,  it  is  predominantly  the  axial 
gradients  of  these  terms  that  impact  the  momentum  and 
energy  balances.  The  radial  momentum  equation  is 
further  influenced  by  the  correlation  terms  arising  from 
the  radial  source  terms  (/.e,,  of  Eqs.  19  and  20). 

The  magnitudes  of  the  correlation  terms  and  their 
axial  gradients  depend  on  the  direction  over  which  the 
governing  equations  are  averaged.  That  is,  the 
correlation  terms  are  not  invariant  through  a  coordinate 
transformation.  Consider  the  flow  fields  schematically 
depicted  in  Figs.  1-3.  For  brevity,  the  discussion  is 
again  restricted  to  constant  radius;  it  is  conceptually 
simple  to  extend  the  analysis  to  account  for  spanwise 
variations  in  the  flow  field.  The  three  scenarios  each 
depict  a  gas  dynamic  wave  that  separates  a  reference 
state  0  which  is  at  rest  relative  to  the  rotor  from  a  state 
1  in  which  flow  is  either  entering  or  exiting  the  rotor 
passage.  If  gas  dynamic  waves  in  wave  rotor  passages 
were  to  always  propagate  with  orientations  that  were 
normal  to  the  blade  surfaces  (as  depicted  in  Figs.  1  and 
2),  then  the  passage-averaging  might  optimally  be 
applied  along  surfaces  oriented  normal  to  the  rotor 
surfaces  (i.e.,  constant- ^  lines  in  Figs.  1  and  2  and 
constant-  ji  lines  in  Fig.  3).  By  averaging  in  the  intrinsic 
(5,1])  coordinates,  the  magnitude  of  the  correlation 
terms  and  their  streamwise  gradients  would  be 
minimized;  indeed,  an  averaging  carried  out  along 
constant- 5  lines  in  an  uncambered  passage  geometry 
(Figs.  1  and  2)  would  lead  to  no  correlation  terms  at  all. 
This  would  hold  both  for  both  steady  and  unsteady 
flows  because  the  transverse  (t]  )  gradients  would  be 
zero  in  the  uncambered  passages.  The  relative  flow 
turning  of  cambered  passages  (e.g.,  see  Fig.  3) 
establishes  transverse  gradients  that  give  rise  to  finite 
correlation  terms,  even  if  the  passage-averaging  is 
carried  out  in  the  intrinsic  coordinates.  The  streamwise 
gradients  of  these  correlation  terms  would  still  be  zero 
in  regions  where  the  flow  is  locally  steady  (Le., 
d/dt  =  0),  for  example,  fore  and  aft  of  the  expansion  fan 
in  Fig.  3.  But,  unlike  the  uncambered  passage,  the 
streamwise  gradients  of  the  correlation  terms  are  non¬ 
zero  in  the  unsteady  flow  regions  of  passages  with 
surface-profiling-induced,  relative,  flow  turning. 

The  above  discussion  might  seemingly  motivate  an 
averaging  procedure  in  intrinsic  coordinates;  however,  in 
general  the  gas  dynamic  waves  do  not  propagate 
oriented  normal  to  the  rotor  surfaces;  rather,  the  wave 


5 


orientation  is  influenced  by  the  gradual  opening  and 
closing  of  the  rotor  passages,  the  passage  geometry,  and 
wave  refraction  caused  by  temperature  and  velocity 
nonuniformities  {c.g.^  at  contact  discontinuities). 
Therefore,  there  seemed  little  impetus  to  formulate  the 
passage-averaged  equations  in  intrinsic  coordinates.  It 
is  noted  however,  that  having  carried  out  the  averaging 
procedure  at  constant  z-lines,  the  gradients  of  the 
correlation  terms  are  also  non-zero  wherever  there  are 
local  span  wise  or  pitch  wise  nonuniformities,  including 
regions  of  steady  flow  turning. 

In  the  following  sections,  the  influence  of  the 
correlation  terms  is  discussed  by  using  analytical 
solutions  for  idealized  gas  dynamic  wave  propagation  in 
representative  wave  rotor  passage  geometries. 

Shock  Wave  in  a  Staggered,  Uncambered  Passage 

Consider  a  shock  wave  propagating  into  a  gas  at 
rest  relative  to  a  wave  rotor  passage  defined  by  the 
staggered,  uncambered  blades  shown  in  Fig.  1.  The  net 
force  established  by  the  shock  wave  is  a  constant  in 

time,  -  VV  ~ 

net  force  is  due  to  the  difference  between  the  pressure 
forces  on  the  pressure  and  suction  sides  of  the  blades  in 
the  instantaneous  "shadow"  of  the  shock  (0  ^  z  ^1). 
By  inspection  it  is  evident  that  the  correlation  terms  are 
non-zero  in  this  same  region  because  the  shock  wave  is 
oriented  at  an  angle  relative  to  the  averaging  surface. 
The  flow  is  locally  steady  (a/3r=0)  outside  of  this 
region  and  the  correlation  terms  are  zero.  The  extent  of 
the  shadow  of  the  shock  relative  to  the  chord  length  is 
sinp/a.  By  considering  that  typical  wave  rotor 
solidities  (o)  are  between  10  and  20  and  that  inflow 
and  outflow  blade  angles  are  as  high  as  |  P  |  ^  45* ,  it  is 
evident  that  the  shock  wave  could  induce  correlation 
terms  over  a  significant  (e.g.,  10%)  portion  of  the  rotor 
passage. 

The  magnitude  of  the  correlation  term  for  two 
arbitrary  variables  /  and  g  in  the  vicinity  of  the  shock  is 
given  by 

3(pS</s'>)  _ 

dz  sinP  cosP  (32) 

The  correlation  terms  are  locally  related  to  the  tangential 
(pitchwise)  differences  ( A^)  in  the  flow  field  variables 
at  the  averaging  surface  which  are  established  by  the 
shock  wave.  The  tangential  difference  terms  in  Eq.  32 
can  be  obtained  by  using  the  following  shock 

relationships:  6(pw  )=w  6(p)  and  8(pw  +/?)  = 

(  SH  C 

w  5(  p  w  ) .  The  relative  magnitudes  of  the  correlation 
ss  ^ 


terms  can  then  be  attained  by  using 

8^(pS<w^w'>)  ^  ^0  ^ 

9(pS<W  >)  ^  ^SH 

z  I  \-y— - ^ - 

(33) 

The  relative  shock  speed,  w  ,  and  the  density  ratio, 

p  /p  ,  are  functions  of  the  shock  strength,  The 

relative  magnitudes  of  the  correlation  terms  vary  with 
shock  strength  and  axial  position  (because  <w^>  varies 
axially).  Evidently,  the  gradients  of  the  correlation 
terms  are  on  the  same  order  of  magnitude  as  the  local 
gradient  of  the  inviscid  flux  in  the  vicinity  of  the  shock 
wave. 


Expansion  from  a  Staggered,  Uncambered  Passage 

Consider  an  expansion  wave  propagating  into  a  gas 
at  rest  relative  to  a  wave  rotor  passage  defined  by 
staggered,  uncambered  blades  shown  in  Fig.  2.  By 
neglecting  end  effects,  the  expansion  wave  propagates 
oriented  normal  to  the  blades  as  shown  and  the  flow 
field  is  uniform  along  the  isobars  shown.  The  simple 
wave  spreads  linearly  with  time;  however,  the  loading 
established  by  the  expansion  fan  is  a  constant  in  time 
(F  L  r  (l-r  /r  ) (1 -€)sinP/0),  like  the  shock 

0  0  i?  7*  H  T 

wave  considered  earlier. 

The  correlation  terms  for  the  uncambered  passage 
are  zero  fore  and  aft  of  the  expansion  fan  where  the 
flow  is  locally  steady  (/.e.,  0/3f  =  O).  The  correlation 
terms  are  non-zero  in  the  regions  of  unsteadiness 
established  by  the  expansion  fan  0  ^  z  ^  1  of  Fig. 
2).  Their  magnitudes  again  depend  on  the  tangential 
differences  at  the  averaging  plane  instantaneously 
established  by  the  expansion  wave  (cf.  Eq.  32).  Unlike 
the  shock  wave,  the  tangential  differences  decrease  with 
time  in  the  case  of  the  expansion  fan  and  must  be 
obtained  from  the  time-dependent  flow  field  solution. 

The  streamwise  velocity  within  the  expansion  fan 
(see  Fig.  2)  is 


w(1,0 

a 

0 


(34) 


where  t  -L  fa  and  where  ^  is  the  position  of  the 
0  t'  0  0 

expansion  fan  at  t  =  0.  The  intrinsic  and  cylindrical 
coordinates  are  related  by 


^  ^0  t  ,  z't  sinp  .  ®  . 

“r'T/  ■  o  ■  6-6^35) 

7  0  0  0  P  S(d5) 
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where  X  jX  =  1  +  - 1]  is  the  ratio  of  the 

wave  speeds  fore  ani  aft  of  the  expansion  fan  and 
where  the  normalized  axial  position  (z  = 
{z-z  )/(z  -z  .  ))  is  defined  in  Fig.  2.  Aft  of  the 
expansion  fan,  the  local  velocity  is 

w  la  - - (X  IX  -1)  which  depends  on  the 

exp^sion  f^  pressure  strength  ( e  =pjp^)^  The  relative 
velocity  in  front  of  the  expansion  fan  is  zero.  The  local 
speed  of  sound  is  related  to  the  local  streamwise 
velocity  component  through 


oce.o  .  v-i"'/?-'* 


1  + 


(36) 


and  the  local  static  pressure,  PlP^^  is  isentropically 
related  to  the  local  speed  of  sound. 

The  relative  magnitude  of  the  correlation  term 
gradients  (e.g.,  d  (pS<wW>)ld  (pS<w^>)  were 
calculated  numerically  by  usmg^  the  ffow  fiel&  described 
above.  The  expansion  fan  spreads  with  time;  therefore, 
the  magnitude  of  the  local  pressure  gradient,  and  hence 
the  flow  field  nonuniformity,  decreases  with  time.  The 
relative  magnitude  of  the  gradient  of  a  correlation  term 
at  a  fixed  normalized  position  (O^z  ^1)  within  the 
expansion  fan  is  a  function  of  the  ratio  of  specific  heats 
(Y),  the  expansion  fan  strength  (e),  the  stagger  angle 
(sinP),  the  passage  solidity  (o),  and  time  (a^r/Lp.  At 
early  times,  the  relative  magnitudes  of  the  correlation 
terms  are  of  order  unity  at  the  ends  of  the  expansion  fan 
and  approach  the  significant  levels  of  the  shock  wave 
considered  in  the  previous  section.  In  general,  the 
relative  magnitudes  of  the  axial  gradients  of  the 

correlation  terms  (3  (p5<w  w^>)/d  (PiS<w^>))  inside 

the  expansion  fan  decrease  with  time  (-  (a^f/Lp  ,  see 

-0  8 

Fig.  4),  decrease  with  solidity  ('-a  '  ,  see  Fig.  5),  and 

increase  with  stagger  angle  (-siii  P  ,  see  Fig.  6).  At 
times  greater  than  the  (normalized)  blade  passing  period 
(a/(rQa)),  the  relative  magnitudes  of  the  axial 
gradients  of  the  correlation  terms  are  small  (0  0.1%) 
except  at  the  head  of  the  expansion  fan  where  the  local 
relative  velocity  is  nearly  zero.  It  is  therefore  expected 
to  be  a  good  approximation  to  neglect  the  correlation 
terms  introduced  by  averaging  procedure  in  a  flow  field 
established  by  a  simple  expansion  fan,  except  at  early 
times  when  the  simple  model  problem  posed  here  is 
invalidated  by  the  passage  gradual  opening  and  closing 
end)  effects. 


Expansion  from  Circular  Arc,  Unstaggered  Passage 
Consider  an  expansion  wave  propagating  into  a  gas 
at  rest  relative  to  a  rotor  passage  defined  by  two 
concentric  circular  arcs  as  shown  in  Fig.  3.  The 
propagation  of  the  expansion  wave  is  a  complicated, 


multi-dimensional  process  which  is  beyond  the  scope  of 
the  present  work  to  describe;  however,  it  is  interesting 
to  consider  the  scenario  depicted  schematically  in  Fig. 
3  in  the  light  of  the  earlier  discussion  of  the  expansion 
fan  in  the  uncambered  passage. 

The  correlation  terms  are  zero  in  the  quiescent  flow 
(state  0)  fore  of  the  expansion  fan.  The  propagating 
expansion  fan  establishes  pitchwise  and  spanwise 
gradients,  and  hence  non-zero  correlation  terms.  The 
magnitudes  of  the  correlation  terms  will  depend  upon 
the  local  acceleration  which  has  two  components  in  the 
case  of  the  cambered  passage:  a.)  the  streamwise,  |i- 
gradient  of  pressure  which  depends  on  the  expansion  fan 
strength  (e)  and  which  is  expected  to  decay  strongly 
with  time  (a  tIL  );  and  b.)  the  transverse,  /^-gradient 
which  is  influenced  by  the  local  flow  turning.  In  the 
circular  arc  geometry,  it  is  expected  that  following  the 
flow  transient  induced  by  the  expansion  fan,  the  radial 
component  of  velocity,  will  soon  relax  to  zero. 
When  >v  =  0 ,  the  flow  behaves  locally  like  a  free  vortex 

R 

flow  in  which  dPIdR  =  pw^lR  and  w  -  IIR.  In  such 

regions  (e,g.,  state  1,  aft  of  the  expansion  fan),  it  is 
expected  that  the  relative  magnitudes  of  the  gradients  of 
the  correlation  terms  will  increase  with  both  the  passage 
radius-ratio  (R  IR  )  and  radius-to-arc-length  ratio  (or 
normalized  radius  of  curvature,  R^IL^).  If  the  circular 
arc  blades  were  staggered,  the  stagger  angle  would 
likely  influence  the  correlation  terms  as  well  through  its 
influence  on  the  instantaneous  "shadow"  of  the  wave. 
Note  that  although  the  flow  is  locally  steady  aft  of  the 
expansion  fan,  because  of  the  averaging  surfaces 
chosen  in  this  work,  and  in  contrast  to  the  staggered, 
uncambered  passage  case,  the  axial  gradients  of  the 
correlation  terms  are  non-zero  in  state  1 :  the  correlation 
terms  impact  the  passage-averaged  momentum  and 
energy  balances  in  the  portions  of  the  passage  in  which 
steady  flow  turning  takes  place. 

Question  of  Need  to  Model  Correlation  Terms 

The  question  of  whether  or  not  it  is  imperative  to 
include  the  correlation  terms  to  attain  the  accuracy 
required  for  engineering  problems  was  not  answered  by 
the  present  study.  Although  beyond  the  scope  of  the 
present  paper,  future  work  should  assess  the  importance 
of  the  correlation  terms  in  the  following  stepwise 
fashion:  first,  while  neglecting  the  correlation  terms,  the 
inviscid  flow  in  a  passage  defined  by  two  concentric 
circular  arcs  should  be  computed  by  using  the  "one¬ 
dimensional"  model  provided  in  the  present  work;  the 
obtained  flow  solution  should  be  compared  with  results 
from  a  multi-dimensional  flow  solver  (e.g..  Ref.  6)  and 
the  accuracy  of  the  one-dimensional  solution  assessed. 
If  the  instantaneous  solutions  of  the  two  solution 
methods  agree,  then  for  the  purpose  of  wave  rotor 
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design  and  analysis  studies,  the  correlation  terms  can  be 
neglected.  If  however  the  results  are  significantly 
different,  then  the  correlation  terms  might  in  fact  be 
relatively  important.  In  this  case,  a  generalized  model 
for  the  correlation  terms  should  be  developed  and 
validated  by  using  results  from  multi-dimensional 
simulations. 

Summary 

The  governing  equations  for  an  inviscid  flow  in 
cylindrical  coordinates  fixed  to  the  spinning  rotor  were 
passage.-averaged  by  integrating  from  hub  to  tip  and 
from  blade  to  blade.  The  averaging  was  motivated  by 
the  need  for  a  reduced  set  of  governing  equations 
appropriate  for  the  parametric  analyses  of  wave  rotor 
design  studies.  Application  of  the  passage-averaged 
equations  requires  closure  models  for  blade,  hub,  and 
tip-shroud  force  integrals  and  for  correlation  terms.  In 
one  of  the  force-integral  models  considered,  the  three 
momentum  equations  were  collapsed  into  a  single, 
effectively  axial,  momentum  equation  by  relating  the 
density-averaged  tangential  and  radial  velocity 
components  to  the  density-averaged  axial  velocity 
component,  under  the  constraint  that  the  flow  travels 
along  a  path  inferred  from  the  rotor  geometry.  This 
approach  yielded  an  essentially  one-dimensional 
formulation  suitable  for  a  CFD-based  analysis  tool  for 
wave  rotors  with  generalized  rotor  passage  geometries. 
It  remains  for  future  work  to  implement  and  to  validate 
this  formulation.  The  validation  can  be  earned  out  by 
comparing  the  predictions  of  the  implemented  model 
with  results  from  a  multi-dimensional  solver. 

Correlation  terms  arise  in  the  application  of  the 
passage-averaged  equations  when  the  density-average  of 
products  of  the  flow  field  variables  are  expanded  into 
products  of  the  density-averaged  variables.  Hie 
magnitudes  of  the  correlation  terms  and  their  gradients 
are  influenced  by  the  orientation  of  the  averaging 
surface.  In  the  present  work,  the  correlation  terms  are 
non-zero  if  there  are  spanwise  or  pitchwise  flow  field 
variations  caused  by  either  local  entropy  gradients  or 
flow  acceleration.  Acceleration  can  be  established  by 
a.)  relative  unsteadiness,  b.)  relative  flow  turning 
induced  by  blade,  hub,  and  tip-shroud  profiling,  and  c.) 
rotation.  The  influences  of  relative  unsteadiness  on  the 
correlation  terms  were  analyzed  by  considering  the 
idealized  problems  of  shock  and  expansion  waves 
propagating  in  staggered,  uncambered  passages.  The 
relative  magnitudes  of  the  axial  gradients  of  the 
correlation  terms  can  be  on  the  order  of  the  axial 
gradients  of  the  inviscid  fluxes  in  the  portions  of  the 
rotor  passage  influenced  by  the  gas  dynamic  waves. 
Relative  flow  turning  induced  by  the  profiling  of  the 
rotor  passage  surfaces  also  establishes  pitchwise  and 
spanwise  gradients,  and  hence  correlation  terms. 
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Figure  !•  Schematic  diagram  of  a  shock  wave 
propagating  into  gas  at  rest  relative  to  a  wave  rotor 
passage  defined  by  staggered,  uncambered  blades 
(blade-to-blade  view). 
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Figure  2.  Schematic  diagram  of  an  expansion  fan 
propagating  into  gas  at  rest  relative  to  a  wave  rotor 
passage  defined  by  staggered,  uncambered  blades 
(blade-to-blade  view). 


Figure  3.  Schematic  diagram  of  expansion  fans 
propagating  into  gas  at  rest  relative  to  wave  rotor 
passages  defined  by  unstaggered,  circular  arc 
(cambered)  blades  (blade-to-blade  view). 


Figure  4.  Relative  magnitude  of  the  gradient  of  the 
correlation  terms  of  the  axial  momentum  equation,  at 
nondimensional  axial  position  (z  )  in  expansion  fan, 
as  a  function  of  nondimensional  time  (f/t^)- 


Figure  5.  Relative  magnitude  of  the  gradient  of  the 
correlation  terms  of  the  axial  momentum  equation,  at 
nondimensional  axial  position  (z  )  in  expansion  fan, 
as  a  function  of  solidity  ( o ). 


sin^p 


Figure  6.  Relative  magnitude  of  the  gradient  of  the 
correlation  terms  of  the  axial  momentum  equation,  at 
non-dimensional  axial  position  (z  *)  in  expansion  fan, 
as  a  function  of  blade  angle  ( P ). 
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